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We reason in support of the universality of quantum spectral fluctuations in chaotic systems, 
starting from the pioneering work of Sieber and Richter who expressed the spectral form factor in 
terms of pairs of periodic orbits with self-crossings and avoided crossings. Dropping the restriction 
to uniformly hyperbolic dynamics, we show that for general hyperbolic two-freedom systems with 
time-reversal invariance the spectral form factor is faithful to random-matrix theory, up to quadratic 
~ order in time. We relate the action difference within the contributing pairs of orbits to properties 

■ of stable and unstable manifolds. In studying the effects of conjugate points, we show that almost 

self-retracing orbit loops do not contribute to the form factor. Our findings are substantiated by 
£S) . numerical evidence for the concrete example of two billiard systems. 
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I. INTRODUCTION 



One of the fundamental questions of quantum chaos is why, in the semiclassical limit, almost all classically hy- 
i— y perbolic systems display universal spectral fluctuations, only depending on their symmetries. This universality was 
initially conjectured by Bohigas, Giannoni and Schmit [1] and is by now supported by overwhelming experimental 
and numerical evidence [2,3]. Examples for experimental tests range from nuclear physics and atomic and molecular 
spectroscopy to classical microwave billiards studied in the limit of large wavenumbers. Even though random-matrix 
theory provides a phenomenological description of these universal features, a derivation from first principles is still 
'—h lacking. Also in other areas of research, such as mesoscopic quantum transport, the reasons for the amazing success 
i " . I of random-matrix theory are only beginning to emerge [4]. To tackle this long-lasting challenge, several approaches 
have been suggested such as parametric level dynamics [2] or an extension of field theoretical methods used in the 
CN ' theory of disordered systems [5,6]. 

In the present paper, following an ansatz pioneered in [7-9], we relate quantum spectral statistics to the properties 
of classical periodic orbits. We treat the spectral form factor (i.e. the Fourier transform of the spectral two-point 
correlation function). According to random-matrix theory, it has the following form for systems belonging to the 
orthogonal universality class (i.e. hyperbolic systems with no symmetries except a time-reversal symmetry whose 
time-reversal operator squares to unity) 
m ■ , 

K(t) = 2t -t ln(l + 2t) = 2r - 2t 2 + 2t 3 . . . . (I) 

Here, r = ^ is the time measured in units of the Heisenberg time Th(E) = 2-Khp(E) (where p(E) is the average 
level density), and we are only considering the range < t < |. Using Gutzwiller's trace formula for the level density 
1* [10], the form factor can be expressed as a double sum over periodic orbits 7, 7' 

X 



K ^ = jr- ( E ex P ( fa ~ s -f ) " Km, - /v)f J s It lx2_izl j \ , (2) 

H \t,7' / E.T 

where A 7 is the classical stability amplitude, T 7 the period and /z 7 the Maslov index of the periodic orbit 7; the 
brackets (. . -)e,t denote averaging over the energy and over a small time window. Note that as the semiclassical 
limit is taken for fixed r, the period of the contributing orbits tends to infinity proportionally to Th- The crucial 
point is that the summand connected to each pair of orbits has a phase given by their action difference divided by h. 
Thus in the semiclassical limit, the phases will be randomly distributed and most terms will interfere destructively A 
contribution to the form factor can only arise from pairs of orbits whose action difference is of the order of Planck's 
constant. Thus, there is a deep relation between correlations among the actions of classical periodic orbits and 
correlations in quantum spectra [8] . 

First success in this direction was reported by Berry [7] , who derived the leading term 2r in the series expansion of 
K(t) from pairs of orbits which are either identical or related by time reversal (diagonal approximation). Starting with 
results by Argaman et al. [8] , the search was on to identify further families of orbit pairs with similar action, expected 
to give rise to higher-order contributions to the form factor (see [11] and references therein). A breakthrough was 
recently achieved by Sieber and Richter, who proposed orbit pairs in which one of the orbits contains a self-crossing 
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in configuration-space with a small angle e [9]. Its partner narrowly avoids that crossing, approximately following one 
loop of the first orbit, and following the other loop in the time-reversed sense (cp. Fig. 1). The action difference 
between the two is quadratic in e and thus can be arbitrarily small. Since it is required that the time reversal of 
a classical orbit loop is again a classical orbit loop, such pairs of orbits can only exist in time-reversal invariant 
systems. Using these pairs, Sieber and Richter derived the leading off-diagonal contribution — 2t 2 to the form factor 
for a uniformly hyperbolic billiard (a billiard where all orbits have the same Lyapunov exponent A), the so-called 
Hadamard-Gutzwiller model, i.e. geodesic motion on a tesselated surface of constant negative curvature of genus 2. 
They summed over the contributions of all these pairs using two ingredients, the action difference between the two 
partners and the density of crossing angles. The r 2 -term arises due to a correction to the latter of next-to-leading 
order in the orbit period. It stems from the fact that in the Hadamard-Gutzwiller model, a loop with a small crossing 
angle must have a minimal traversal time £ m i n (e) = — j logce. In their derivation these authors made use of several 
system-specific niceties of the Hadamard-Gutzwiller model, which the form factor, being universal, cannot depend 
on. First steps towards an extension to other systems have been taken for quantum graphs in [12,13], where both the 
t 2 - and the r 3 -contribution to the spectral form factor were shown to originate from similar orbit pairs. 

In the present paper, extending results we first presented in [14,15], we go beyond these idealized models and 
derive the leading off-diagonal contribution to the form factor for general two-dimensional hyperbolic systems with a 
Hamiltonian of the form i?(Q,P) = ^ + V(Q). The main novel ideas needed to establish this universality are: 

(i) The relation between the partner orbits can be formulated elegantly in terms of the invariant manifolds, which 
also determine the action difference within the orbit pair. 

(ii) The Maslov indices of the partner orbits can be shown to coincide. 

(iii) In general systems, a logarithmic correction to the angle distribution arises involving the Lyapunov exponent of 
the system (as defined below). 

(iv) In systems with conjugate points, the one-to-one correspondence between crossings and orbit pairs is broken. 
There are crossings related to almost self-retracing loops without an associated partner orbit, and "braids" of crossings 
with a common partner. We show how to overcome these problems and reveal the distribution of crossings for which 
a partner orbit does exist as universal even in the presence of conjugate points. 

(v) The universal contribution to the form factor follows from a relation between the invariant manifolds and the 
Lyapunov exponent. 

Our findings are substantiated by numerical results for two billiard systems. 

Note that the use of crossings in configuration space constitutes no conceptual problem e.g. concerning canonical 
invariance, as what we are using is in fact the geometry of the invariant manifolds. As pointed out in [16], it appears 
natural to work in configuration space since it is singled out by the conventional time-reversal operator. Interestingly, 
the results of [16] imply that our considerations immediately carry over to systems with non-conventional time-reversal 
invariance, as long as they can be canonically transformed to a Hamiltonian of the above structure. An alternative 
approach avoiding a projection to configuration space will be presented in [17,18]. 

This paper is organized as follows. We first determine the action difference between the two partner orbits (Sect. 
2) and show that their Maslov indices coincide (Sect. 3). In Section 4 we investigate the statistics of crossings angles, 
phase-space locations and loop times in systems without conjugate points. In Section 5 we generalize our findings to 
systems with conjugate points and clarify the relation between crossings and orbit pairs in such systems. Finally, we 
will derive the r 2 -contribution to the form factor in Section 6. 



II. ACTION DIFFERENCE 

A long periodic orbit typically has a huge number of self-crossings in configuration space. For a given small-angle 
self-crossing, we will show that a related orbit with similar action exists, provided the two orbit loops separated by 
the crossing are long 1 . This partner orbit is obtained from the initial orbit by time reversal of one loop and a local 
deformation close to the crossing. Further away, the deviation between the initial orbit (respectively its time-reversed) 
and the partner decays exponentially. As proposed in [16], the relation between these two orbits can be expressed 
elegantly in terms of stable and unstable manifolds. The main advantage of this method is that the action difference 
within the orbit pair can be shown to depend only on the crossing angle and the local behavior of the invariant 
manifolds at the phase-space location of the crossing. The complicated dependence on stability matrices of the two 
loops derived in [9] is thus drastically simplified in the limit of long loops. 



1 Certain subtle issues concerning the existence and uniqueness of that partner will be dealt with in Section 6. 
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Recall that two phase-space points lie on each other's stable (unstable) manifold if trajectories starting from these 
points come infinitely close for t — > oo (t — > — oo). For two- freedom systems, the local behavior of these invariant 
manifolds at a phase-space point X can be characterized by just one number, called "direction" or "curvature" of 
the given manifold [20,21]. This direction is defined as the ratio between the momentum and configuration-space 
components p and q of an infinitesimal phase-space deflection along the stable (unstable) manifold 



dp 
dq 



(X); 



(3) 



here, the deflection is assumed orthogonal to the trajectory passing through X. Note that time reversal T exchanges 
stable and unstable manifolds and changes the sign of their directions such that B S (K) = — _B„(TX). 

Now, our condition of a local deformation can be formulated in terms of the invariant manifolds. Let Xi and X2 
denote the phase-space points corresponding to the two traversals of the crossing. On the partner orbit with time- 
reversed right loop, we single out a phase-space point X' 1; located where the perpendicular from the crossing hits that 
orbit (see Fig. 1). X' : must approximately lie on the stable manifold of Xi (and vice versa), as trajectories starting 
there approach for a long time as t — > 00. Conversely, X^ must lie on the unstable manifold of the time reversal of X 2 . 
To determine X^, it is convenient to work in a Poincare section S defined by the above perpendicular. The coordinates 
of a phase-space point X = (Q, P) in that section will be denoted by x = (q,p). Since — xi = (q[ — qi,p[ — pi) 
must be stable and x[ — Tx 2 = (q[ — qi,p'\ + Pi) unstable, we obtain the following system of linear equations in q[, 
p[ involving the directions of the invariant manifolds 



pi-pi =B a (Xi)(gi - qi ) 
p' 1 +P2=B u {K' 1 )(q[-q 2 ). 

Furthermore, we know that Xi and X 2 coincide in configuration space, i.e. 
orthogonal momenta of Xi and TX 2 is related to the crossing angle by p\ -f 
(4) can be written as 



|Pi|e 



B u (Xi)-B a (Xi) 



(4) 

also qi = g 2) and that the difference in 
P2 = — |Pi|e. Using this, the solution to 



(5) 



Pi follows trivially. 




FIG. 1. Sketch of a Sieber-Richter pair in configuration space. Phase-space points are denoted by arrows starting at the 
corresponding configuration-space location, and momentum is indicated by the direction of the arrow. Depicted are the 
phase-space points of the two traversals of the crossing and two phase-space points of the partner orbit as defined in the text. 
The Poincare sections orthogonal to X' x and X' 2 and two points half-way through the loops divide the orbit into four parts a, 
b, c, and d. 



One can now expand the action of the orbit containing the crossing around that of its partner avoiding it. As the 
action of a periodic orbit is stationary, the first-order term vanishes. We restrict ourselves to the quadratic order in 
Qi — <?i,2- We first consider the "upper" side of the crossing, i.e. the orbit parts a and b separated by S and two 
points half-way through the two loops as shown in Figure 1. The second derivatives of their action can be related 
to the stable and unstable directions. For example, S a (Q' 1 ,Q' L ) generates the motion from the configuration-space 

point Qj on S to the point half-way through the left loop, thus P[ = —§^r and therefore = ~^r- These 

derivatives are taken for constant Q' L , and can in the limit of long loops be approximated by derivatives along the 
stable manifold. We thus have 

(Xi) = -B a (Xi), (6) 
and analogously 



d 2 S a _ dp[ 
dq' 2 dq[ 
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d 2 S b = d P [ 

dq[ 2 dq[ 



(X' 1 ) = B U (X[). (7) 



The same reasoning can be repeated for the "lower" side of the crossing in a slightly different Poincare section 
orthogonal to the "lower" part of the orbit in a phase-space point X' 2 (see Fig. 1). 

A Taylor expansion now shows that avoiding a crossing with angle e located at X reduces the action by 

p2 2 p2 2 

2(B U (X[) - B a (Xi)) + 2(B U (X^) - B a (X£)) ' w 

For small angles, we can neglect the difference of the stable resp. unstable directions at Xi, X^, TX 2 , and TX' 2 
(which we will denote collectively by X), as long as B s and B u are sufficiently smooth close to the crossing location. 
This condition is usually fulfilled unless e.g. immediately after the crossing, the different branches of one loop enclose 
a singularity of the flow. Apart from these exceptional cases, (8) simplifies to 

A ^> = g .<x)-k<xr (9 » 

This may also be negative, meaning that the action can also be increased by avoiding the crossing. Using that 
B S (X) = — B U (TX) one easily sees that this result is invariant under time reversal of X, as it should be, since we 
must obtain the same action difference if we insert for X the other, almost time-reversed traversal of the crossing. 
Assuringly, the action difference in the Hadamard-Gutzwiller model follows from (9) as a special case. In that model 
we have B U (X) = —B S (X) = mX (where m is mass) for all phase-space points X, and thus AS* = in accordance 
with [9]. 

III. MASLOV INDEX 

The contribution of a pair of orbits to the form factor is determined not only by the difference of their actions, 
but also by the difference of their Maslov indices. We want to show that for Sieber-Richter pairs of orbits, at least 
in the limit of small angles the latter vanishes. Note that this was trivial in the Hadamard-Gutzwiller model, where 
all Maslov indices are zero. We work in the framework of a beautiful geometric interpretation of the Maslov index 
of a periodic orbit due to Creagh et al. [2,22], again in a (in the present context two-dimensional) Poincare section 
orthogonal to the orbit. Here, the invariant manifolds locally have the form of lines through the origin which rotate 
around the origin as we move along the orbit. The Maslov index now equals the net number of clockwise half- 
rotations of the stable or, equivalcntly, the unstable manifold (i.e. the difference of the numbers of clockwise and 
counter-clockwise half- rotations) . 

For our argument, we also define the Maslov index of a finite non-periodic orbit such as an orbit loop in Sieber- 
Richter theory. It is the sum of the net rotation angles of the stable and the unstable manifold around the origin, 
divided by 2ir. Even though not canonically invariant, non-integer, and depending on the units chosen in the Poincare 
section, this definition is very useful because it makes the Maslov index of a loop invariant under time reversal. 
Obviously, time reversal leaves the absolute value of a rotation angle invariant. The same is true for the sense of 
rotation, since time reversal inverts the motion on the Poincare section in direction (turning a clockwise rotation into 
a counter-clockwise one and vice versa), but also changes the sign of the momentum (turning the sense of rotation 
back to the original one) 2 . An alternative analytic proof can be found in [23]. In addition, time reversal exchanges 
the stable and unstable manifolds, which also cannot affect the sum of their rotation angles. Reassuringly, this newly 
defined Maslov index coincides with the usual one in case the orbit happens to be periodic. Moreover, it is additive 
for subsequent orbit pieces and smooth under small deformations of the orbit as long as the invariant manifolds are 
smooth. 

We can now show that the Maslov indices of the two partner orbits coincide. They can be expressed as sums over 
the Maslov indices of the two loops. Due to additivity and time-reversal invariance of the loop Maslov indices, formal 
time reversal of one loop leaves the Maslov index of the orbit invariant. In the limit of small crossing angles, any 
subsequent local deformation as described in Section 2 can at most lead to a small change of the Maslov indices of 
the loops. Since, however, the Maslov index of a periodic orbit is an integer quantity, the Maslov indices of the two 
partners have to coincide. 



2 Note that time reversal also inverts the directions of both coordinate axes in the Poincare section, but this has no impact on 
the sense of rotation. 
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IV. CROSSINGS IN SYSTEMS WITHOUT CONJUGATE POINTS 



We have seen that in non-uniformly hyperbolic systems the action difference within a Sieber-Richter pair of orbits 
depends both on the angle and on the phase-space location of the crossing. Anticipating that like in the Hadamard- 
Gutzwillcr model, also the traversal times of the loops play a crucial role, we investigate the density p(K,e,t\T) of 
loops with crossing angle e, time t, and initial phase-space point X (the latter being one of the two traversal points 
of the crossing) in a periodic orbit of period T. This density is normalized in such a way that integration over all 
possible values of X, e, and t yields the average number of loops, and thus twice the average number of crossings, 
in a periodic orbit of period close to T. To make our exposition more clear, in the present Section we will still limit 
ourselves to systems without conjugate points; our findings will be generalized to systems with conjugate points in 
Section 5. 

From the ergodicity of the flow we can, along the lines of [9] , deduce the following approximation 

2P 2 

p crg (X,e,t|T) = -^Tsine, (10) 

where |0| is the volume of the energy shell. The proportionality to sine reflects the fact that orthogonal parts of the 
orbit intersect more frequently than almost parallel or antiparallel ones. 

However, like in the Hadamard-Gutzwiller model there is a correction to that ergodic prediction because loop times 
below a certain angle and location dependent minimum are impossible. We will first give a general argument why this 
minimal loop time is generic, and later substantiate our point by numerical evidence for billiards and discuss certain 
system-specific limitations. In systems without conjugate points, two trajectories starting from the same point in 
configuration space with a small opening angle cannot gather in another point as long as the phase-space separation 
between them can be approximated as a linear function of the initial separation. For an orbit loop to close it is thus 
necessary that, during half the loop time the separation between one of its two branches and the time-reversed of 
the other branch grows far enough to make that approximation invalid. Apart from exceptional cases to be discussed 
later, this typically means that their separation must reach a classical phase-space scale of the order of some finite 
fraction of the maximal phase-space separation. If the branches enclose a small crossing angle e, this requires the loop 
to be long. Let ||<5X|| denote the phase-space separation between one branch of the loop and the time-reversed of the 
other one in an arbitrary norm |.| at the location of the crossing (thus ||5X|| oc e) and 5X(|) the separation after 
half the loop time. As long as the linearized approximation is applicable both are for large t related by 
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(11) 



where A is the Lyapunov exponent of the system. This asymptotic behavior is due to the fact that long loops ergodically 
explore the whole energy shell. The Lyapunov exponent governing the asymptotic fate of a small deviation 5X (with 
non-vanishing unstable component) at a phase-space point X coincides for almost all X, since it can be expressed as 
the average of the so-called local stretching rate over an infinite trajectory starting at X [20]. Thus, in case of an 
ergodic flow it coincides with a phase-space average almost everywhere (except e.g. on periodic orbits). Note that due 
to Pesin's theorem, A also coincides with the Kolmogorov-Sinai entropy [20,24]. In uniformly hyperbolic systems such 
averaging is trivial since the local stretching rate is constant and therefore even the Lyapunov exponents of periodic 
orbits coincide. 

Demanding that the deviation ||<5X(|)|| reaches the limit for a breakdown of the linear approximation just after 
half the loop time, we obtain a minimal loop time of 

t min (X,e) = -|lnc(X)e. (12) 

Here, c(X) is constant with respect to the angle but may depend on the phase-space point X immediately preceding 
the loop. However, the exact value of c(X) will turn out to be irrelevant for the form factor. The second loop starts 
approximately at TX, thus its time must fulfill T — t > i m ; n (TX, e). Incorporating these minimal loop times in a 
straight-forward way, we obtain the following density of loops 

p(X, e, t\T) = Pcrg (X, e, t\T)Q(t - t min (X, e))6(T - t - i min (TX, e)). (13) 

Integration over < t < T yields the density of crossing angles and phase-space locations 

2P 2 

P(X, e\T) = ^^T(T - t min (X, e) - t min (TX, e)) sin e. (14) 
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We want to clearly point out the scope of our approximation. First, the minimal loop time derived above gives 
the threshold after which loops have a chance to close, but does not predict the exact behavior of the shortest loops. 
Thus in the immediate vicinity of the minimal loop time, system specific structures appear which are not described 
by (13) and (14). Second, we will see that singularities can let the linear approximation for the separation between 
the branches of a loop break down even before a typical phase-space scale is reached (as also discussed in [17]). Thus, 
singularities give rise to exceptional, system-specific crossings which may ignore the minimal loop time [14]. 

For our numerics, we are also interested in the statistics of self-crossings of non-periodic orbits. In full analogy to 
the above considerations for periodic orbits, we can show that the density of loops with crossing angle e, time t, and 
initial phase-space point X in non-periodic orbits with traversal time T reads 

2P 2 

p n P(X,e,i|T) = -j^p(T-t)e(t-t min (X,e))sine. (15) 
Integration yields the density of crossing angles and phase-space locations 3 

P n P(X,e|T) = J^(T-£ min (X,e)) 2 sine. (16) 



A. Example: The desymmetrized diamond billiard 

We will now present numerical evidence for this crossing distribution for a special billiard system, the desymmetrized 
diamond billiard (see Fig. 2). It can be regarded as the empty space between four overlapping disks [25] a la Sinai, cut 
into eight equal pieces. The distance between the disks is chosen as one, and we choose their radius r = 0.541 so that 
its interior angles become § , \ and J . The desymmetrized diamond billiard is non- uniformly hyperbolic and belongs 
to the class of semi-dispersing billiards (i.e. it is surrounded by a boundary which consists of locally concave and 
straight segments). Therefore, it is free of conjugate points. It has a circumference C = 0.671 and an area A — 0.0157. 
Santalo's formula [26] gives its mean free path as I = ^ = 0.0735. By averaging over the Lyapunov exponents of 
random non-periodic trajectories, we numerically obtain the Lyapunov exponent of the system as A = 4.31. 



FIG. 2. The desymmetrized diamond billiard (with an example for a singularity-related crossing). 

A few words are in order about our numerical technique. The crossing statistics is determined by averaging over 
2 x 10 7 non-periodic trajectories of length L = 10 (i.e. very long orbits compared to the typical length scales I and 
j) with random initial conditions. For the delicate statistics of angles < yi, even 10 9 trajectories were considered. 
Note that for billiards, it is useful to work in dimensionless coordinates with mass and velocity equal to one; then 
the traversal time T of each orbit equals L, and the same holds true for the the traversal times t and lengths I of the 
loops. 

We now turn to our results. For small angles, the basic idea of a minimal loop time depending logarithmitically 
on the crossing angle can be verified at a glance at Figure 3a. Here we have plotted the density of angles and loop 
times (up to t = 3) of crossings occurring anywhere on the energy shell 4 , and for convenience divided out the term 
sin e arising from the ergodic approximation. Sufficiently far above that logarithmic threshold, marked by a dashed 
line in Figure 3a, this density is almost uniform. In agreement with (15), it decays linearly towards larger loop time. 



3 Note that in contrast to periodic orbits, for each crossing there is only one loop, hence it is only included once in this 
distribution. 

4 In contrast to the following results, for this density plot it was sufficient to only take into account 5 x 10 7 trajectories. 
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Note that according to (12), the minimal loop time weakly depends on the location on the energy shell due to the 
factor c(X). Thus, the minimum in Figure 3a is slightly smeared out. Below the dashed line, the density of crossings 
diminishes fast before vanishing completely inside the white region. 




FIG. 3. Density plot of the combined distribution of loop times t and crossing angles e in the desymmetrized diamond 
billiard: a) for e < 0.057T, b) for all angles. The density is normalized such that multiplied by sin ededt, it gives the number of 
crossings in the respective intervals of the crossing angle e and the loop time t inside one orbit of traversal time T = 10; the 
resulting scale is shown in c). For small angles, we observe a threshold logarithmic in e, as indicated by a dashed line. 

In addition to this expected behavior, a rich variety of system-specific structures are seen. As announced, the 
crossing density shows system-specific inhomogeneities in the immediate vicinity of the minimal loop time. In addition, 
below the minimal loop time, we observe exceptional crossings related to singularities of the flow, most importantly 
the tangential singularity. Namely, the linear approximation for the separation between the two branches of a loop is 
already violated if one branch reflects e.g. at the circular part of the boundary, and the other branch narrowly avoids 
the circle. In this case, the separation between the branches of the loop need not reach a typical phase space scale 
for the loop to close. An example for such exceptional loops ignoring the minimal loop time is shown in Figure 2. 
However, our numerical results indicate that their effect on the crossing distribution is minute. 

At larger angles, which in the semiclassical limit give no contribution to the form factor, further structures are seen 
(cp. Fig. 3b), such as 

(i) Discrete lines at e = n corresponding to periodic orbits, which are nothing but loops with a "crossing angle" n. 
They are deformed and broadened when going to smaller angles. Thus, most loops close to the minimal time are 
obtained by deformation of the shortest periodic orbits. 

(ii) Four families of loops which start at zero loop time, and as shown in [14] are close to corners. One of these families 
of loops involves only reflections at the two straight-line segments of the boundary and has a crossing angle of J for 
all its members. 

Note that most of these structures were absent in the idealized example of the Hadamard-Gutzwiller model. In 
that model, the analog of Figure 3 just consists of dispersionless logarithmic curves corresponding each to the family 
of loops obtained by deformation of one periodic orbit [9] . 

Further evidence for the minimal loop time is shown in Figure 4, this time also taking into account the phase-space 
dependence of the crossing statistics. We restrict ourselves to crossings which have angles smaller than a maximal 
angle e max and take place inside a bin whose area is | of the total energy shell, and consider the density of their loop 
times t. This density must have a gap for small t whose width is given by the minimal loop time and thus increases 
logarithmitically when e max is reduced. Incidentally, this density of loop times can be shown to decay linearly for large 
t in the case of non-periodic orbits. Our numerical results confirm these predictions and again reveal system-specific 
structures close to the minimal loop time. 

The angle density P np (X, e\T) (averaged over the whole energy shell and over bins with sizes of | and of 
the energy shell) agrees with our predictions as well (cp. Fig. 5). For sufficiently small angles, the distribution is 
sinusoidal. Upon division by sine, in accordance with (16) a logarithmic correction due to the minimal loop time 
becomes visible. The relative weight of that correction (seen as the slope in the logarithmic plot in Fig. 5b) coincides 
for all three samples. We thus see that indeed the crossing distribution in a billiard is homogeneous on the energy 
shell apart from system-specific oscillations for large e and apart from the factor c(X). Based on the relative weight 
of the logarithm, fitting yields a minimal loop time i m i n (X, e) = — In c(X)e where Afi t = 4.31 agrees perfectly with 
the Lyapunov exponent A. 
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loop time t loop time t 

FIG. 4. Statistics of loop times in non-periodic trajectories of traversal time T = 10 in the desymmetrized diamond billiard. 
Depicted is the density of crossings with loop time t and crossing angle < e max in a fraction of | of the energy shell for a) 
£max = 7T, b) e max = 0.l7r, c) e max = 0.0l7r, d) e max = 0.0017T. The result is normalized such that multiplied by dt it gives 
the number of such crossings in the loop time interval (t, t + dt) inside one orbit of traversal time T; the graphs are based on 
averages over 2 x 10 7 orbits for large angles and 10 9 orbits for small angles, respectively. Due to the minimal loop time, there 
is a gap for small t which g rOWS elS Cm ax 

is reduced. 
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FIG. 5. a) Average of P np (X, e\T) over the whole energy shell of the desymmetrized diamond billiard and fractions of | and 
thereof (divided by for normalization), b) The same quantities divided by ^irp in a logarithmic plot. The logarithmic 

correction due to the minimal loop time becomes visible. 



V. CROSSINGS IN SYSTEMS WITH CONJUGATE POINTS 



We now want to generalize our treatment to systems with conjugate points. In such systems, trajectories fanning 
out from the same point in configuration space with a small opening angle can focus again in a second point, then in a 
third, etc. (see Fig. 6). All these points are called mutually conjugate. Wc will sec that conjugate points destroy the 
one-to-one relation between crossings and orbit pairs, because there are (i) crossings without a partner orbit avoiding 
the crossing and (ii) families ("braids") of crossings with a common partner. Both effects have a direct analogy to 
quantum graphs, for which the leading off-diagonal contribution to the form factor was derived in [12] and even the 
third order of the expansion was obtained [13]. The general picture emerges that a periodic orbit has one partner for 
each two almost time-reversed orbit stretches dividing the orbit into sufficiently long loops. 
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FIG. 6. Example for a family of mutually conjugate points in the cardioid billiard. 

Even though our findings are general, we find it instructive to also discuss their meaning in the special case of 
systems with symbolic dynamics. Here, each periodic orbit is unambiguously defined by a string of symbols. Each 
symbol in the alphabet denotes one partition of a certain Poincare section, and the symbol sequence of a periodic 
orbit is composed of the symbols corresponding to the partitions it traverses. 

In particular, we consider the example of the cardioid billiard, which belongs to the family of focusing billiards 
(i.e. billiards surrounded by a locally convex boundary). The cardioid has been intensively studied in the literature, 
see [27-32] and references therein. It is hyperbolic, and the fidelity of both symmetry-reduced spectra to random- 
matrix theory was demonstrated numerically in [28] . It has symbolic dynamics with two symbols effectively denoting 
straight-line segments of the orbit [29-31]. The initial points of each segment and the cusp divide the boundary into 
two parts. A symbol A is assigned to the segment if, seen from the initial point, the final one lies in the part on the 
clockwise side, and a symbol B if it lies on the counter-clockwise side. We note that time reversal inverts the ordering 
of symbols and interchanges A and B. For a more detailed account of the following results, we refer the reader to [14]. 

A. Crossings without partner 

In systems with conjugate points, loops below the minimal loop time derived above are possible, because the two 
branches of a loop can meet while the linear approximation for their separation is still applicable. However, we will 
see that in this case the orbit has no partner avoiding the crossing. 

Examples for these loops are shown in Figure 7. Here, a loop starts from the crossing, is reflected with an almost 
right angle and nearly retraces itself before crossing itself with a small angle. The locations of the crossing and the 
reflection are almost conjugate to each other, since the two traversals of the crossing limit a fan of trajectories with 
a small opening angle which gathers again at the reflection point. 




FIG. 7. Examples for crossings conjugate to a reflection in the cardioid billiard. 

These loops are so close to being time-reversal invariant that the "partner" formally determined by time reversal 
of one such loop coincides with the orbit itself. To show this, we use the notation introduced in Section 2 and let 
L and R denote the stability matrices of the left and right loop separated by the crossing (cp. Fig. 1). For small 
crossing angles, the partner with time-reversed right loop has to fulfill the following linear system of equations (given 
in different form in [9]) for its traversals x^, x.' 2 of a Poincare section orthogonal to the first traversal of the crossing 5 



5 We might as well consider the slightly more complicated Poincare sections we had to choose in Section 2 for technical reasons. 
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x 2 - x 2 = L(x[ - Xl) 
xi -Tx 2 = i? r (x 2 -Txi), 



(17) 



where i? r = Ti? X T is the stability matrix of the time-reversed right loop. Note that in contrast to the reasoning in 
Section 2, here we do not require the orbit loops to be long. With <5xi = Xj — xi, Sx 2 = x 2 — x 2 , 5x = Tx 2 — xi = 
(0, |Pi|e), (17) simplifies to 



This system of equations in Sxi and 5x 2 is valid up to corrections quadratic in <Jx, and obviously has exactly one 
solution. By going to higher orders, Spehner indeed rigorously showed that there is exactly one "partner orbit" [17]. 
It is, among all periodic orbits, unambiguously singled out by fulfilling (18) in linear order. This "partner orbit" 
does however coincide with the initial orbit, if the right loop is smaller than the minimal loop time derived above. In 
this case, the separation between the two branches of the loop can be treated in a linear approximation. During the 
traversal time of the right loop x 2 is carried into xi, and Txi is carried into Tx 2 . Since the separation between Tx 2 
and Xi can be treated in a linear approximation using the stability matrix of the right loop or its time-reversed, we 
obtain (up to quadratic order in <5x) 



Comparing this to the second equation in (18), we see that (18) has the trivial solution <Jxi = <5x 2 = 0. Thus, the 
"partner" with time-reversed right loop coincides with the initial orbit. Conversely, if the left loop were shorter than 
the minimal loop time, that "partner" would coincide with the time reversal of the initial orbit. 

In both cases, as pairs of identical or mutually time-reversed orbits are already included in the diagonal approx- 
imation, such crossings related to almost self-retracing loops give no off-diagonal contribution to the form factor. 
We conclude that in systems with conjugate points, the minimal loop time obtains a new meaning: A partner orbit 
avoiding a given crossing exists only if both loops separated by that crossing exceed the minimal loop time. In the 
sequel, we will refer to crossings with an associated partner orbit as "relevant", and to the others as "irrelevant" 
crossings. The crossing statistics derived above for systems without conjugate points immediately carries over to 
relevant crossings in systems with conjugate points. We refer to the Appendix for a discussion of the statistics of 
irrelevant crossings. 

For systems with symbolic dynamics, the irrelevance of almost self-retracing orbit loops can be shown even easier. 
Each crossing divides the symbol sequence of the orbit in two parts L and R corresponding to the two loops 6 . The 
partner can be determined by reverting the symbol sequence of one loop in time. An orbit with symbol sequence LR 
thus has a partner with symbol sequence LR , where the superscript r denotes time reversal. The two are identical 
up to time reversal if one of the symbol sequences L or R is time- reversal invariant. This provides a rigorous criterion 
to decide whether for a given crossing, there is an orbit avoiding the crossing or not. 



We want to present numerical evidence that the distribution of relevant crossings in the cardioid indeed conforms 
to (15) and (16). We looked for self-crossings in 8 x 10 6 non-periodic orbits of length L = 250 (which is long compared 
to the mean free path 1.851 [30] and the Lyapunov length 2.83); for crossings with angles < y^, even 4 x 10 8 orbits 
were considered. In order to exclude irrelevant crossings from our statistics, we used the symbolic-dynamics criterion 
derived above. 

Again, we have to deal with the effect of singularities on the crossing statistics. In the cardioid billiard, there are 
system-specific crossings related to the cusp, which will be discussed in more detail in the Appendix. Since for the 
corresponding orbits, even the applicability of the Gutzwiller trace formula is questionable, we find it necessary to 
also distinguish in our numerics between generic and cusp-related crossings. Like irrelevant crossings, the latter are 
excluded from our statistics using a symbolic-dynamics criterion derived in the Appendix. 



For simplicity, we assume here that each loop can unambiguously assigned a symbol sequence. Minor technical difficulties 
arise for symbols which effectively denote the crossing segments and not parts of a loop; they are dealt with in [14]. 



5x 2 = LSxi 
Sxi -5x = R T (Sx 2 +TSx). 



(18) 



R T TSx = R T (x 2 - Txi) = xi - Tx 2 = -Sx. 



(19) 



B. Numerical results 
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Our numerical results for generic relevant crossings display a striking similarity to the case of the desymmetrized 
diamond billiard, thereby clearly supporting our point that the observed effects are universal. Again, the combined 
density of crossing angles and loop times 7 reveals a minimal loop time depending logarithmitically on the angle, and 
some system-specific inhomogeneities close to that minimum (cp. Fig. 8). Furthermore, the distribution of loop times 
corresponding to crossings in a given energy-shell bin and with a maximal angle e max shows a gap for small angles 
(see Fig. 9). Finally, the density P np (X, e\T) of generic relevant crossings shown in Figure 10 (averaged over the 
whole energy shell and over fractions of | and ^ thereof) agrees well with (16). Based on the relative weight of the 
logarithm we obtain a fitting value for the Lyapunov exponent of Afi t = 0.352 coinciding up to a minute numerical 
error with A = 0.353 8 . 




FIG. 8. Density plot of the combined density of loop times t and crossing angles e in the cardioid billiard: a) for e < 0.057T, 
b) for all angles. Normalization as in Figure 3; the resulting scale is shown in c). For small angles, we observe a threshold 
logarithmic in e, as indicated by a dashed line. 



a) 




b) o. 




C) o.a 




loop time t 



d) ,c 




loop lime 1 



FIG. 9. Statistics of loop times in non-periodic trajectories of traversal time T = 250 in the cardioid billiard. Depicted is 
the density of crossings with loop time t and crossing angle < e max in a fraction of | of the energy shell for a) e max = -k, b) 
£max = 0.1-7T, c) e max = 0.017T, d) e m ax = 0.0017T. Normalization as in Figure 4; the graphs are based on averages over 8 x 10 6 
orbits for large and 4 x 10 orbits for small angles. As in Figure 4, due to the minimal loop time there is a gap for small t 
which grows as e max is reduced. 



7 In contrast to the following results, for this distribution it was still sufficient to include 3.5 x 10 7 trajectories. 
s This value follows from the results of [30] for the Kolmogorov-Sinai entropy of the billiard map ft map = IX = 0.653 and could 
be reproduced by averaging over the Lyapunov exponents of non-periodic orbits. 
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crossing angle £ crossing angle £ 

FIG. 10. a) Average of P np (X, e|T) over the whole energy shell of the cardioid billiard and fractions of | and <rj thereof 
(divided by ^r^p for normalization), b) The same quantities divided by in a logarithmic plot. The logarithmic correction 
due to the minimal loop time becomes visible. 

C. "Braids" of crossings with common partner 

In systems with conjugate points, two almost mutually time-reversed stretches of an orbit cross several times and 
thus have a whole family ("braid") of small-angle crossings located close to mutually conjugate points. This is because 
two trajectories starting at the traversals of a given crossing will meet again in points conjugate to it, thereby forming 
new crossings. An example for such a braid of crossings in the cardioid is shown in Figure 11a. 




FIG. 11. a) Example for braid of crossings in mutually conjugate points in the cardioid billiard. We depict the central part 
of such a braid occurring in a periodic orbit, b) Crossings in the corresponding partner orbit. The number of crossings differs 
exactly by one. 

Interestingly, no matter which of these crossings one tries to avoid, one always obtains the same partner orbit. To 
understand this, consider two crossings close to mutually conjugate points (cp. Fig. 12). We introduce Poincare 
sections orthogonal to one traversal Xi of the first and one traversal Xi of the second crossing. Furthermore, let L, 
R (L, R) be the stability matrices of the loops separated by the first (second) crossing. If M is the stability matrix 
describing the motion from Xi to Xi, we have (in the notation introduced above) 

L = M T LM 

R = M- 1 R(M T )- 1 

<5x = M-^x. (20) 

The partner corresponding to the first crossing intersects the Poincare section orthogonal to that crossing at x^ = 
Xi + Sxi , and a linear approximation shows that it intersects the Poincare section orthogonal to the second crossing 
at x£ = Xi + <5xj with 

<5xi = M _1 <5xi 

<5x 2 = M r <5x 2 . (21) 
It is easy to show that these Siti also fulfill the system of equations determining the partner which corresponds to 
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the second crossing (i.e. the "tilded" version of (18)) 9 . Consequently, for each braid of crossings close to mutually 
conjugate points and thus for any two approximately time-reversed stretches of an orbit there is just one partner 
orbit. 

R 



FIG. 12. Sketch of a Sieber-Richter pair with crossings in mutually conjugate points. The arrows denote phase-space points 
on the Poincare sections defined in the text. L, R, L, R, and M are stability matrices as described in the text. 

Again, we would like to point out how these findings translate to symbolic dynamics. The existence of two almost 
time-reversed stretches of an orbit is naturally reflected in its symbol sequence, because the stretches will have mutually 
time-reversed symbol sequences Z and TJ . Thus, small-angle crossings appear in orbits with symbol sequences of the 
form IZrZ r (where I and r are arbitrary symbol sequences) and occur between the stretches belonging to Z and Z r 
[19]. In case of conjugate points, there are several such crossings between these stretches. Any of them divides Z into 
two parts Z| and Z r , and the whole symbol sequence into parts L = Z\ T \Z\ and R = Z r rZ r T corresponding to the loops. 
The partner obtained by time reversal of the right loop has the symbol sequence 

LR r = Z| r IZ|(Z r rZ r r ) r = Z| r IZ|Z r r r Z r r = IZr r Z r , (22) 

regardless of which crossing was chosen. Here, we used that time reversal of a symbol sequence inverts the ordering 
of its subsequences and reverts the subsequences in time, and that symbol sequences related by cyclic permutation 
are equivalent. Thus, we see again that the orbit has only one partner for the whole braid. 

This partner, too, contains a braid of crossings, as demonstrated for the example of the cardioid billiard in Figure 
lib. For both orbits within the pair, the crossing angles increase approximately exponentially towards the edges of the 
braids, where the two orbit stretches deviate more and more from being mutually time-reversed [14]. The crossings 
of the partner orbit are slightly shifted compared to those in the original orbit (most strongly close to the center), 
and the numbers of crossings in both orbits differ by one, the orbit with larger action containing one crossing more. 

This observation, which will turn out to be crucial for the derivation of the form factor, will be proved in the sequel 
for general hyperbolic Hamiltonians of the form H(Q, P) = + V(Q). Note that it also trivially applies to systems 
without conjugate points, where the partner with larger action contains one crossing and the other none. Our proof 
relies on an argument of winding numbers. We follow one of the two almost time-reversed stretches of the orbit and 
study, in a Poincare section orthogonal to the orbit, three quantities, the stable and unstable manifolds (which locally 
have the form of straight lines through the origin) and the (small!) phase-space vector <5x pointing to the time reversal 
of the other orbit stretch (cp. Fig. 13). As we move along the orbit, these lines and vectors rotate around the origin. 
Every time 5x rotates through the p-axis, the orbit has a crossing. Note that for kinetic-plus-potential Hamiltonians 
the p-axis may only be traversed in clockwise direction, because because we have q — ^ > in the upper and < in 
the lower half plane. 

For the orbit with larger action, our formula for the action difference demands that whenever a crossing occurs we 
have -B U (X) — £> S (X) > 0, i.e. the unstable manifold has a higher slope in the Poincare section than the stable one. 
Thus for the partner with larger action, Sx is located between the stable manifold (on the counter-clockwise side) 
and the unstable manifold (on the clockwise side), and for the partner with smaller action, the inverse is true. The 
motion of 5x is given as a superposition of the rotation of invariant manifolds and a motion from the stable towards 
the unstable manifold (since the deviation between the orbit stretches first shrinks and then increases). We thus see 
that for the orbit with larger action, <5x performs one more clockwise half- rotation around the orbit. Consequently, it 
crosses the p-axis once more and therefore the orbit with larger action contains one more crossing 10 . 



9 A similar result holds in a phase-space based treatment, as will be shown in [17]. 

10 This result holds even true in case of hard-wall reflections. Here, Sx does not rotate continuously but jumps, and q changes 
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FIG. 13. Poincare section orthogonal to one ol two almost time-reversed orbit stretches with stable and unstable manifold 
and the vector 5x pointing to the time reversal of the other stretch. Here, <5x traverses the p-axis, i.e. the stretches cross in 
configuration space. For both Sieber-Richter partners, the asymptotic motion of 5x with respect to the invariant manifolds is 
indicated by arrows. 



VI. LEADING OFF-DIAGONAL CONTRIBUTION TO THE FORM FACTOR 




Once the crossing statistics and the action difference are known, we can now determine the contribution to the 
form factor arising from Sieber's and Richtcr's family of orbit pairs. It follows from (2) that each pair of orbits of 
period T = tTh gives a contribution of 2 A 2 cos ' to the form factor (where we neglect the difference in amplitude 
and period between the two partners). However, we have to make sure that in spite of the one-to-one correspondence 
between crossings and orbit pairs being lost in general systems, each pair of orbits is counted only once. The key to 
the solution is to formally assign to each crossing a contribution of 2 A 2 cos ^psign(i? u (X) — £? S (X)). Since the partner 
with larger action and thus B U (X.) — B S (X.) > contains one more crossing that the partner with B„(X) — B S (X.) < 0, 
these formal contributions of all crossings in both partners add up to the correct value for the pair. We can thus sum 
over all Sieber-Richter pairs by summing over all orbits with period T = tTh and integrating over the locations X of 
the crossings on the energy shell ft and over their angles e 

V^(T-T 7 ) / d 3 Xsign(B u (X)-B s (X)) [" de P(X,e|T) cos A5( *' e) \ . (23) 
~ Jn Jo n I E.T 

Note that each crossing is counted twice, since it is traversed at two mutually time-reversed phase-space points. As 
there are also two partner orbits with time-reversed left or right loop, no additional factor appears here. The integral 
over the angle can be performed in a way similar to [9] using (9), (12) and (14) and yields in the semiclassical limit 

LB fX) — B 

— t 2mA|n| ^ n * s vame om y depends on the asymptotic action difference and angle distribution for e — > 0. 
Interestingly, the contribution of the ergodic approximation of the angle distribution vanishes, and the result is only 
due to its logarithmic correction originating from the minimal loop time. Applying the sum rule of Hannay and Ozorio 
de Almeida [33] (£ 7 A 2 S(T - T 7 )) T = T, we are led to 

- -^^r 1 - < 24 » 

where (...) denotes an average over the energy shell. 

Now, ergodic theory comes into play, relating the invariant manifolds to the Lyapunov exponent. The (positive) 
Lyapunov exponent of a hyperbolic two-freedom system can be expressed as the energy-shell average of the so-called 
local stretching rate [20] 



X (X) = tr 



d 2 H d 2 H 



dQdP d 2 P 



C„(X) 



(25) 



sign without traversal of the p-axis. However, as the two partners undergo the same number of hard-wall reflections, this does 
not affect the difference in numbers of crossings. 
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Here, C M (X) is a matrix which relates the momentum and configuration-space components of unstable deviations by 
dP = C„(X)c£Q. For systems with Hamiltonian iJ(Q,P) = |^ + V(Q), this local stretching rate is proportional to 
£?„(X). Evaluating the trace in coordinates orthogonal and parallel to the orbit, we see that x(X) = ^trC u (X) = 
^B U (~K). Thus, we infer that (B u ) = mX. Proofs for the special case of semi-dispersing billiards can also be found 
in [34]. Due to B S (TX.) = —B U (X.), we also have (B s ) = —(B u ) = —m\. From this, we immediately obtain the 
universal leading off-diagonal contribution to the spectral form factor 

K 2 (t) = -2r 2 . (26) 

Acknowledgements: I am grateful to Peter Braun, Fritz Haake, Stefan Heusler, and Dominique Spehner for 
close cooperation. Furthermore, I want to thank Arnd Backer, Sven Gnutzmann, Gerhard Knieper, Klaus Richter, 
Holger Schanz, Martin Sieber, Lionel Sittler, Uzy Smilansky, Marko Turek, Carlos Viviescas, and Wen-ge Wang for 
enlightening discussions. 



VII. APPENDIX: MORE ABOUT CROSSINGS IN THE CARDIOID BILLIARD 



Crossings in reflection-conjugate points 



It is instructive to see how the statistics of angles, locations and loop times of irrelevant crossings in "reflection- 
conjugate" points (i.e. points conjugate to the location of a self- retracing reflection) compares to the statistics of 
relevant crossings. We again consider the cardioid, where the distinction between these classes of crossings can be 
made with the help of symbolic dynamics. 

Numerically, we observe that the density of angles belonging to irrelevant crossings reaches a constant value for 
e — > (sec Fig. 14a), while for relevant crossings the corresponding density was shown to be proportional to sine up 
to a logarithmic correction. Thus, we see that in fact the vast majority of small-angle crossings are irrelevant. 

In configuration space, irrelevant crossings are concentrated in the immediate vicinity of certain curves, the loci 
of reflection-conjugate points (see 14b). One can construct these curves by traveling along the boundary, starting 
trajectories in a direction perpendicular to the boundary, and determining points conjugate to these starting points 
using simple geometric optics. One such curve (sometimes divided into several parts) exists for the first, second, etc. 
conjugate points met after the starting point. 

For each of these loci, the length (and thus time) of the trajectory from the boundary point to the crossings 
(corresponding to half of the length of the self-retracing loop) has several local maxima as shown in Figure 14c. A 
short calculation shows that these maxima are responsible for peaks in the loop time distribution which decay sharply 
towards larger and smoothly towards shorter time [14] (see Fig. 14d). 
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FIG. 14. Statistics of irrelevant crossings with almost self-retracing loops in the cardioid billiard, a) Density of angles e 
corresponding to irrelevant crossings in non-periodic orbits of traversal time T = 250 (divided by ). b) Distribution in 

configuration space. A random sample of irrelevant crossings found numerically is marked by crosses. They cluster close to 
the loci of reflection-conjugate points. The locus of the first (second) conjugate point met after each self-retracing reflection is 
marked by a solid (dashed) curve, c) Traversal times of the self-retracing loops, parametrized by the arc length of the reflection 
point on the boundary. Again, the different curves refer to the first, second, etc. conjugate point met after the reflection point, 
d) Density of loop times for loops with crossing angle < yjy . 

A. Cusp-related crossings 

In this Appendix, we want to show that in the cardioid billiard, the minimal loop time is ignored by a system-specific 
class of crossings related to the singularity at the cusp. We will argue that for the evaluation of the form factor, these 
crossings have to be disregarded as well. 

An example for these crossings is shown in Figure (15). For one of the two loops separated by the crossing, the two 
branches stay very near until being reflected on opposite sides of the cusp. Thus, the linear approximation for the 
separation between these branches breaks down before their separation reaches a typical phase-space scale. Then, in 
many cases, the loop closes by going around the boundary like in a whispering gallery The length (and thus time) 
of the shortest such loops is approximately given by the circumference of the billiard. In fact, they can be seen as 
deformations of finite orbits starting and ending at the cusp as introduced in [29,30]. In spite of their short time, the 
crossing angle can be arbitrarily small. So for these special loops, there is no minimal time depending logarithmitically 
on the angle. 




FIG. 15. An example for a cusp-related loop and its symbol sequence. The loop here has a symbol sequence of the form 
SATBAS r (with S = B and T = B 6 ). Similar loops following the patterns SBTABS T , SABTAS r , and SBATBS T are obtained 
by time reversal and reflection at the symmetry line. 

However, they do give rise to pairs of orbits with similar action, as their symbol sequences are not time-reversal 
invariant. Nevertheless, we have to disregard these crossings for a semiclassical evaluation of the form factor, as even 
the validity of the Gutzwillcr formula for orbits coming so close to the cusp is questionable. In general, each time 
a periodic orbit comes closer to the boundary of a billiard than some length scale of the order of Planck length, 
so-called penumbra corrections to the trace formula play an important role, as shown in [35] for dispersing billiards. 
The breakdown of the trace formula for families of orbits approaching the cusp of the cardioid is discussed in [29,32]. 

Due to these problems, we have to distinguish in our numerics between generic and cusp-related crossings. This 
distinction is done based on symbolic dynamics (compare Fig. 15). The loop parts preceding the cusp are almost 
time-reversed with respect to each other and thus have symbol sequences S and S r . In contrast, the part following 
the cusp is almost symmetric with respect to the symmetry line of the cardioid, i. e. its symbol sequence T has to be 
invariant under inverting the ordering of symbols [30]. Usually, it undergoes a series of subsequent either clockwise 
or counter-clockwise reflections at the boundary, i.e. T just consists of several identical symbols A or B. In addition, 
for the symbols between S, T and S r one has to take into account the effect of one branch being reflected close to the 
cusp and the other one narrowly avoiding it. These symbols can be read off from the example shown in Figure 15. 
Thus, we see that cusp-related crossings have symbol sequences of the form SATBAS r , SBTABS r , SABTAS r , and 
SBATBS r , where S is arbitrary and T invariant under inverting the ordering of symbols. 



16 



[1] O. Bohigas, M. J. Giannoni, C. Schmit, Phys. Rev. Lett. 52, 1 (1984) 

[2] F. Haake, Quantum Signatures of Chaos, 2nd edition , Springer, Berlin (2000) 

[3] H.-J. Stockmann, Quantum Chaos: an introduction, Cambridge University Press, Cambridge (1999) 
[4] K. Richter, M. Sieber, Phys. Rev. Lett. 89, 206801 (2002) and references therein 
[5] O. Agam, B.L. Altshuler, A. V. Andreev, Phys. Rev. Lett. 75, 4389 (1995) 
[6] R. S. Whitney, I. V. Lerner, R. A. Smith, Waves Random Media 9 (1999), 179 
[7] M. V. Berry, Proc. R. Soc. Lond. A, 400, 299 (1985) 

[8] N. Argaman, F. Dittes, E. Doron, J. Keating, A. Kitaev, M. Sieber and U. Smilansky, Phys. Rev. Lett. 71, 4326 (1993) 
[9] M. Sieber, K. Richter, Physica Scripta T90, 128 (2001); M. Sieber, J. Phys. A 35, 613 (2002) 

[10] M.C. Gutzwiller, J. Math. Phys. 12, 343 (1971) 

[11] U. Smilansky, B. Verdene, J. Phys. A 36, 3525 (2003) 

[12] G. Berkolaiko, H. Schanz, R. S. Whitney, Phys. Rev. Lett. 88, 104101 (2002) 
[13] G. Berkolaiko, H. Schanz, R. S. Whitney, J. Phys. A 36, 8373 (2003) 
[14] S. Miiller, Diplomarbeit, Universitat Essen (2001) 

[15] S. Miiller, talk given at Complex Quantum Systems seminar, Regensburg (2002) 
[16] P. A. Braun, F. Haake, S. Heusler, J. Phys. A: Math. Gen. 35, 1381 (2002) 
[17] D. Spehner, J. Phys. A 36, 7269 (2003) 
[18] M. Turek, K. Richter, J. Phys. A 36, L455 (2003) 

[19] P. A. Braun, S. Heusler, S. Miiller, F. Haake, Eur. Phys. J. B 30, 189 (2003) 

[20] P. Gaspard, Chaos, scattering, and statistical mechanics, Cambridge University Press (1998) 

[21] N. Chernov, R. Markarian, Introduction to the Ergodic Theory of Chaotic Billiards, IMCA, Lima, 

http://premat.fing.edu.uy (2001) 
[22] S. C. Creagh, J. M. Robbins, R. G. Littlejohn, Phys. Rev. A42, 1907 (1990) 
[23] J. A. Foxman, J. M. Robbins, J. Phys. A: Math. Gen. 30 8187 (1997) 

[24] Ya. B. Pesin, Math. USSR Izvestija 10, 1261 (1976); Ya. B. Pesin, Russ. Math. Surv. 32, 55 (1977) 
[25] R. Artuso, J. Tech. Phys. 37, 273 (1996) 

[26] L. A. Santalo, Integral Geometry and Geometric Probability, Addison Wesley, Reading, Mass. (1978) 

[27] M. Robnik, J. Phys. A: Math. Gen. 17, 1049 (1984) 

[28] A. Backer, F. Steiner, P. Stifter, Phys. Rev. E 52, 2463 (1995) 

[29] H. Bruus, N.D. Whelan, Nonlinearity 9 1023 (1996) 

[30] A. Backer, H.R. Dullin, J. Phys. A.: Math. Gen. 30, 1991 (1997) 

[31] A. Backer and N. Chernov, Nonlinearity 11, 79 (1998) 

[32] A. Backer, Dissertation, Universitat Ulm (1998) 

[33] J. H. Hannay, A. M. Ozorio de Almeida, J. Phys. A: Math. Gen. 17, 3429 (1984) 

[34] Ya. G. Sinai, Russ. Math. Surv. 25, 137 (1970); Ya. G. Sinai, Harvard Univ. Preprint (1978); N. Chernov, Funct. Anal. 

Appl. 25, 204 (1991); N. Chernov, J. Stat. Phys. 88, 1 (1991) 
[35] H. Primack, H. Schanz, U. Smilansky, I. Ussishkin, Phys. Rev. Lett. 76, 1615 (1996); same authors, J. Phys. A: Math. 

Gen. 30, 6693 (1997); see also references therein for other semiclassical approaches to diffraction 



17 



